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Introduction 



The decays of pseudoscalar mesons to light leptons are hehcity suppressed in the Standard Model 
(SM) due to the V—A structure of the charged current coupling. In particular, the SM width of 
— >• decays with P = tt, K, D, B (denoted as P12 below) is 

where G_p is the Fermi constant, Mp and Mi are the meson and lepton masses, fp is the 
meson decay constant, and Vgg, is the corresponding CKM matrix element. Although the SM 
predictions for P£2 decay rates are affected by hadronic uncertainties via the decay constant, 
ratios of decay rates of the same parent meson do not depend on fp and can be computed 
very precisely. In particular, the SM prediction for Rk = T{Ke2)/^iKfj_2), inclusive of internal 
bremsstrahlung (IB) radiation, is [1] 

= {Vtj )'(^ ^^'^Q-) = P.47T±0.001) X 10-. 

where 6Rqed = (—3.79 it 0.04)% is the electromagnetic correction. 

Within extensions of the SM involving two Higgs doublets, Rk is sensitive to lepton flavour 
violating effects induced by loop processes with the charged Higgs boson (H^) exchange [2]. 
A recent study [3] has concluded that Rk can be enhanced by 0(1%) within the Minimal 
Super symmetric Standard Model. However, the potential new physics effects are constrained 
by other observables such as Bs — )• fJ'^fJ'~ and i?+ — t- t"*"/^ decay rates [1]. Moreover, Rk is 
sensitive to the neutrino mixing parameters within SM extensions involving a fourth generation 
of quarks and leptons [5] or sterile neutrinos [6] . 

Measurements of Rk have recently been reported by the KLOE [7] and NA62 [8] experiments. 
An improved measurement based on the full dedicated data sample collected by the NA62 
experiment in 2007-2008 and superseding the earlier result [8] is reported here. 

1 Beam and detector 
1.1 Beam line 

The beam line of the earlier NA48/2 experiment [9] was used for the NA62 data taking in 
2007-2008. Either simultaneous or single beams of positive and negative secondary hadrons, 
with central momentum of 74 GeV/c and momentum spread of ±1.4 GeV/c (rms), were derived 
from the primary 400 GeV/c protons extracted from the CERN SPS and impinging on a 40 cm 
long, 0.2 cm diameter beryllium target. The beam momenta were selected by the first two 
magnets in a four dipole achromat and by momentum-defining slits incorporated into a 3.2 m 
thick copper/iron proton beam dump, which also provided the possibility of blocking either of 
the two beams. The beam composition was dominated by pions (vr^), with kaon {K^) fractions 
of about 6%. The and K~ beams entered the decay fiducial volume at angles of ±0.23 mrad 
(±0.30 mrad in the early stage of data taking, about 25% of the total beam flux) with respect to 
the detector axis, so as to compensate for the opposite ±3.58 mrad deflections by the downstream 
spectrometer magnet. These deflections were regularly reversed during the data taking. The 
individual beam particles were not tagged, and their momenta were not measured. The beam 
kaons decayed in a fiducial volume contained in a 114 m long cylindrical vacuum tank. 

The hadron beams were accompanied by an intense flux of stray muons travelling outside the 
beam vacuum pipe. Two 5 m long magnetized iron toroids with small horizontal and vertical 
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apertures centered on the beam line were installed upstream of the decay volume to suppress 
backgrounds associated with these "halo" muons. These toroids, named "muon scrapers" , were 
operated with the same magnetic field polarity chosen to deflect positive halo muons away from 
the beam region, thereby generating a strong charge asymmetry of the muon halo. 

1.2 Detector 

The momenta of charged decay products were measured by a magnetic spectrometer, housed 
in a tank filled with helium at approximately atmospheric pressure, placed downstream of the 
decay volume. The spectrometer comprised four drift chambers (DCHs), each consisting of 8 
planes of sense wires, and a dipole magnet located between the second and the third DCH which 
gave a horizontal transverse momentum kick of 265 MeV/c to charged particles. The measured 
spectrometer momentum resolution was (Tp/p = 0.48% © 0.009% • p, where the momentum p 
is expressed in GeV/c. A counter hodoscope (HOD) consisting of two planes of orthogonal 
plastic scintillator strips producing fast charged particle trigger signals was placed after the 
spectrometer. 

A 127 cm {27Xq) thick liquid krypton (LKr) electromagnetic calorimeter, used for lepton 
identification and as a photon veto detector in the present analysis, was located further down- 
stream. Its 13248 readout cells had a transverse size of 2x2 cm^ each with no longitudinal 
segmentation. The energy resolution was aE/E = 3.2%/\/!E © © 0.42% {E in GeV). The 

spatial resolution for the transverse coordinates x and y of an isolated electromagnetic shower 
was fja; = o-y = 0.42 cln/^/E © 0.06 cm {E in GeV). 

A more detailed description of the detector components used for this measurement can be 
found in Ref. [10]. 

1.3 Trigger logic 

A relatively low beam intensity (corresponding to ~ 10^ kaon decays in the vacuum tank per 
second) was used to enable the operation of a minimum-bias trigger configuration with high 
efficiency, and to minimize the accidental background. The K(,2 trigger condition consisted of 
the coincidence of signals in the two HOD planes (the Qi signal), loose lower and upper limits on 
the DCH hit multiplicity (the 1-track signal) and a LKr energy deposit of at least 10 GeV (the 
EhKr signal). The K^2 trigger condition required a coincidence of the Qi and 1-track signals 
downscaled by a factor of -D = 150. The 1-track condition was not used in the early stage of the 
data taking (about 10% of the total beam flux). Downscaled control triggers were collected to 
monitor the performance of the main trigger signals. 

2 Data samples 

The data used for this measurement were obtained from about 3.5 x 10^ SPS spills (with ~ 10^^ 
protons per spill), collected in 4 months of operation, and correspond to about 2 x lO^'' 
decays in the vacuum tank. The data-taking strategy was optimized to measure the two main 
backgrounds in the i^e2 sample, which are due to the beam halo muons and to K^2 decays with 
a muon (/i^) misidentified as an electron (e^). 

As no kaon tracking is available, beam halo muons are a direct source of background to 
K^2 decays, as well as to decays via — )• e^fgf^ decays in the fiducial region {v^ is used 
to denote either a neutrino or an antineutrino here and below). The muon scrapers installed 
in the beam line were optimized for halo background suppression in the -ftT^ data samples (as 
quantified below), making the decays more favourable for the measurement. To measure 
the muon halo background directly from data, the and K~ data samples were collected 
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alternately by blocking the negative or the positive beam, respectively. Therefore, 65% (8%) of 
the total 2007 beam flux corresponded to {K~) decays collected in single-beam mode. In 
addition to being the signal samples (i.e. providing the K12 data), these data sets are used as 
control samples to measure the muon halo background to the decays of opposite sign kaons (see 
Section [3.3.ip . The remaining 27% of the 2007 beam flux corresponded to decays collected 
with simultaneous beams with the ratio of kaon fluxes of ^(K~^)/^{K~) ~ 2, and cannot be 
used for the halo background subtraction. An additional K~ data sample collected in 2008, 
corresponding to about 4% of the total 2007 beam flux, is used for halo subtraction in the 
sample, but not for the Rk measurement. 

To estimate the K^2 background, the probability to misidentify a muon as an electron due to 
large energy deposition in the LKr calorimeter has been measured. This required the collection of 
a muon sample free from the typical ~ 10~^ electron contamination due to //^ — >• e^Ueiy^ decays 
in flight. To this end, 55% of the kaon flux in 2007 was collected with a transverse horizontal 
lead (Pb) bar installed below the beam pipe between the two HOD planes, approximately 1.2 m 
in front of the LKr calorimeter. The bar was 9.2Xq thick in the beam direction (including an 
iron holder) and shadowed 11 rows of LKr cells (about 10% of the total number of cells). For 
a 50 GeV electron traversing the Pb bar, the probability of depositing over 95% of its initial 
energy in the LKr is ~5 x 10~^, as estimated with a simulation. 

Due to the different acceptance and background conditions, and K~ decays, as well as 
data collected with and without the Pb bar, are analyzed separately. The four resulting inde- 
pendent data samples are denoted as K~^(Ph), K~^{iioPh), K~{Fh) and K~(noPb). The earlier 
analysis [8] is based on the {noFh) data set only, which contains 41% of the reconstructed 
Ke2 candidates and has the lowest background contamination. The present analysis extends to 
the whole data sample and involves several improvements on the estimation of backgrounds and 
systematic uncertainties. 



3 Data analysis 

3.1 Analysis strategy 

The analysis is based on counting the numbers of reconstructed Ke2 and i^^2 candidates col- 
lected simultaneously; therefore it does not rely on an absolute kaon flux measurement. As 
a consequence, several systematic effects cancel to flrst order. Due to the dependence of the 
acceptance and background on the lepton momentum, the Rk measurement is performed in- 
dependently in 10 lepton momentum bins covering a range from 13 to 65 GeV/c (the lowest 
momentum bin spans 7 GeV/c, the others are 5 GeV/c wide). Since the /K~ and Pb/noPb 
samples are treated independently, the analysis is performed separately for 40 statistically in- 
dependent subsamples with partially correlated systematic uncertainties. The ratio Rk in each 
subsample is computed as 

N{Ke2)-N^{Ke2) A{K^2) U X J_ 

D ' N{K^2) - Nb{K^2) ' A{Ke2) ' fe X e(/^e2) " /LKr' 

where N{K£2) are the numbers of selected K£2 candidates {i = e, fi), N-q{Ki2) are the numbers of 
background events, A{K ^2) / A{Ke2) is the ratio of the geometrical acceptances (the acceptance 
correction), fn are the lepton identiflcation efficiencies, e{Ki2) are the trigger efficiencies, /LKr 
is the global efficiency of the LKr readout (affecting the Ke2 reconstruction only), and D = 150 
is the K^2 trigger downscaling factor. 

To evaluate the acceptance correction and the geometrical parts of the acceptances for most 
background processes entering the computation of N-b{K£2), a detailed Monte Carlo (MC) sim- 
ulation based on GeantS [11] is used that includes the time variation of data-taking conditions. 
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The particle identification, trigger and readout efficiencies, as well as the muon halo background, 
are measured directly from data. The determination of each term in the above expression is 
discussed below. 

3.2 Event reconstruction and selection 

Charged particle trajectories and momenta are reconstructed from hits and drift times in the 
spectrometer using a detailed magnetic field map. Fine calibrations of the spectrometer field in- 
tegral and DCH alignment are performed by monitoring the mean reconstructed — )■ Tr^vr+vr" 
invariant mass and the mean reconstructed missing mass in i^^2 decays. Clusters of energy de- 
position in the LKr calorimeter are found by locating the maxima in the digitized pulses from 
individual cells. Shower energies are corrected for energy outside the cluster boundary, energy 
lost in isolated inactive cells (0.8% of the total number) and cluster energy sharing. The en- 
ergy response is calibrated with samples of electrons and positrons from K — t- decays. 
Further details about the reconstruction procedures can be found in Ref. [lOj . 

Most selection criteria are common to both the and Kfj_2 decay modes. The principal 
criteria are the following. 

• Exactly one reconstructed charged particle track (lepton candidate) geometrically con- 
sistent with originating from a kaon decay is required. The geometrical consistency is 
determined by reconstructing the decay vertex as the point of closest approach of the 
lepton candidate track extrapolated upstream (taking into account the measured stray 
magnetic field in the vacuum tank) and the axis of the kaon beam of the corresponding 
charge (determined with fully reconstructed — )• tt^tt^tt" decays). The reconstructed 
closest distance of approach (CDA) of the lepton track to the beam axis is required to 
be less than 3.5 cm, which is crucial for the suppression of the muon halo background, as 
well as backgrounds from decays followed by pion or muon decays in flight. The CDA 
computation also deflnes the decay vertex. 

• The track impact points in the DCHs, HOD and LKr calorimeter must be within the 
corresponding fiducial acceptances, including appropriate separations from detector edges 
and inactive LKr cells. Moreover, it is impossible to efficiently identify electrons traversing 
the Pb bar by energy deposition in the LKr. Therefore, the region of the LKr calorimeter 
shadowed by the bar is excluded from the geometrical acceptance for the K^(Ph) data 
samples, which leads to a reduction of the acceptance by 18% for both and Kf^2 decays 
for these samples. 

• The reconstructed lepton momentum must be in the range 13 to 65 GeV/c. The lower 
limit ensures high efficiency of the Ei^kv energy deposit trigger condition (see Section [L3]) . 
Above the upper limit, the analysis is affected by larger systematic uncertainties due to 
backgrounds, as most backgrounds discussed in Section [331 increase at high track momen- 
tum. 

• No LKr clusters with energy above 2 GeV and within 12 ns of the track time are allowed, 
unless they can be associated to the track via direct energy deposition or bremsstrahlung. 
This requirement provides a photon veto to suppress backgrounds from — t- e^uj, 

and — >• vr^Tr" decays. In the i^^(Pb) data samples, the LKr clusters 
located in the shadow of the bar are not used for the veto condition, as the LKr response 
for photons traversing the bar is difficult to be reproduced by simulation. This reduces 
the photon veto coverage and increases backgrounds from the above decays, as quantified 
in Section [3.3.51 
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• The reconstructed kaon decay vertex must be located within the vacuum decay volume: its 
longitudinal coordinate ^vertex must satisfy the condition ^vertex > -Zmin- Here Zmin depends 
on the reconstructed lepton momentum and is optimized for the suppression of the muon 
halo background, as described in Section r3.3.1l and illustrated for the ir^(noPb) sample 
in Fig. [T^,b. This requirement removes ~ 95% (~ 80%) of the halo background in 
{K~2) samples and decreases the geometrical acceptance for events by about 25%. 

• The residual muon halo background affecting the -fC^ samples (Fig. [TJa) is found to have 
specific geometrical properties and can be strongly reduced by suitable cuts. The points 
(xc, yc) defined by extrapolation of the lepton track to the final collimator plane (a trans- 
verse plane located at the beginning of the vacuum decay volume, z = 0) have a localized 
distribution, as shown for the -fCj[^(noPb) sample in Fig. [TJ:. The optimization of a com- 
mon selection condition for K~2 and K~2 decays (as required to minimize the bias on Rk) 
has been driven by the distribution of the background in the K~2 sample, which is higher 
than the one in the K~2 sample. The regions with highest population of background to 
the K~2 decay indicated in Fig. [lb are rejected, which reduces the halo background in the 
K~2 samples by about 75%, while decreasing the geometrical acceptance for K^2 events by 
about 10%. 

Kinematic identification of the K12 decays is based on the reconstructed squared missing 
mass assuming the track to be an electron or a muon: M'^^^^{t} = {Pk — Pe)'^, where Pk and 
P^ {I = e, fi) are the kaon and lepton 4-momenta, with mi used to assign the components of 
Pi. Pk is taken from the average beam momenta, which are monitored with fully reconstructed 
— )■ TT^TT^TT decays. Fig. [2^ shows the squared missing mass M^-gg(e) evaluated in the 
electron hypothesis for and -fC^2 events as a function of lepton momentum. A selection con- 
dition — < M^jgg(^) < Ml is applied. The limits and Ml vary across lepton momentum 
bins, taking into account the resolution on M^^^^{i), radiative tails and background condi- 
tions. M^ varies from 0.013 (GeV/c^)^ in the central region of the lepton momentum range to 
0.016 (GeV/c^)^ at low and high momenta. Similarly, M| varies from 0.010 to 0.013 {GeV/c^f 
for K^{noPh) samples and from 0.010 to 0.011 (GeV/c^)^ for K^{Ph) samples. The latter 
limits are stricter to compensate for the weaker photon veto. 

Lepton identification is based on the ratio E/p of energy deposition in the LKr to momentum 
measured by the spectrometer, as illustrated in Fig. [2)3. Charged particles are identified as 
electrons if {E/p)mm < E/p < 1.1, where (-E/p)min = 0.95 for p > 25 GeV/c and {E/p)ram = 0.9 
otherwise. The relaxed condition at low lepton momentum is possible because backgrounds in 
the Ke2 sample due to /i^ and vr^ misidentification are rejected kinematically in that momentum 
range (as seen for the Kf^2 background in Fig. [2^). For {E/p)^[^ = 0.95, this criterion leads 
to an electron identification efficiency fe > 99% and a probability of misidentifying a muon as 
an electron of ~ 4 x 10~^. Charged particles with E/p < 0.85 are classified as muons; the 
corresponding muon identification inefficiency is negligible (1 — ~ 3 x 10~^). 

3.3 Backgrounds 

3.3.1 Muon halo background in the samples 

The rate and kinematical distribution of the muon-halo-induced backgrounds are qualitatively 
reproduced by a dedicated simulation. However, the precision of this simulation is limited by the 
uncertainties on the fringe fields of the beam line magnets. Therefore, these backgrounds have 
been measured directly by reconstructing the {^12) decay candidates from control samples 
of positive (negative) tracks collected with the positive (negative) beam blocked, as described 
in Section [2j The control samples used for background subtraction in the four independent K£2 
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Figure 1: Properties of the muon halo background to i^^(noPb) decays measured using control 
data samples as explained in Section [3. 3. 11 Distributions of the reconstructed (a) e+ and (b) e~ 
candidate momentum versus the reconstructed longitudinal coordinate of the kaon decay vertex. 
The momentum-dependent Zmin parameter optimized for muon halo suppression is indicated by 
solid lines: events to the left of the lines are rejected. The residual background in the K~2 sample 
is a factor of 5 higher than in the sample, (c) Distribution of the e~ candidate crossing 
points (xc, Uc) at the final collimator plane {z = 0) for the halo background in the -ftr^(noPb) 
sample surviving the Zmin cut. The regions enclosed within solid lines are excluded to suppress 
the background. Halo background distributions for the -fr^(Pb) samples are similar. 
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Figure 2: (a) Reconstructed squared missing mass in the electron mass hypothesis M^;gg(e) as a 
function of lepton momentum for K(^2 and K^2 decays (data). The wrong mass assignment for 
the K^2 decays leads to the momentum-dependence of M^jgg(e). (b) E/p spectra of electrons 
and muons (data) measured from — )• -K^e^u and K^2 decays. The part of the muon spectrum 
above E/p = 0.95 is for the muons traversing the Pb bar. The electron identification criterion 
applied for p > 25 GeV/c is indicated with arrows. 
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signal samples are mutually exclusive, leading to independent statistical uncertainties. These 
samples are normalized to have the same numbers of muon events with —0.3 (GeV/c^)^ < 
-^miss(/") ^ —0.1 (GeV/c^)^ as the data (such events are not compatible with a kaon decay and 
can only arise from a halo muon). 

As a cross-check, the backgrounds in the Kf>2 samples have also been evaluated with a 
hybrid MC technique, using the measured spatial, angular and momentum distributions of the 
halo muons and simulating their decays (/i^ — >• e^i^e^/x) assuming unpolarised muons, as the 
polarization is unknown. The results agree with those obtained from the direct measurements. 

The non-zero probability of reconstructing a Ke2 candidate due to the decay of an opposite 
sign kaon enhances the halo background estimates for the K(,2 samples. The effect is more 
pronounced for the K~ samples because the ratio of the to K~ beam fluxes is ~ 2. A 
decay must result in at least three charged daughter particles to produce an opposite sign 
particle. Contributions from — t- tt^j^I^u, — )• vr^vr^, and — t- l^ve^e~ [12] decays 
(where i = e, fi, and vr^ — t- e~^e~j denotes the Dalitz decay) have been identified and subtracted 
using MC simulations. The corresponding correction to the final result is negligible: ARk /Rk ~ 
10-4. 

In addition to being charge asymmetric (as explained in Section 11.11 and shown in Fig. [T^, 
b), the muon halo background is left-right asymmetric (as illustrated in Fig. \l}p for the K~2 
sample), and therefore depends on the polarity of the spectrometer magnetic field. To reduce 
the combined statistical uncertainty of the halo background estimates, almost equal samples 
of signal and control data were taken with each spectrometer polarity. The residual polarity 
imbalance was corrected by assigning weights to the control samples. 

As control data have been collected mostly without the Pb bar, parts of these Er^(noPb) 
samples are used to estimate background in the K^(Ph) signal samples by reducing the geomet- 
rical acceptance as described in Section 13. 2i This minimizes the overall statistical uncertainty. 
The background to signal ratio in the region covered by the Pb bar is lower than in the region 
outside the Pb bar due its geometrical localization. Therefore the backgrounds in the K^(Ph) 
samples are higher than those in the i^^(noPb) samples. 

The Ki2 selection procedure and geometrical acceptance are time-dependent due to the 
variations of the beam geometry and the presence of temporarily masked groups of LKr cells. 
This implies a possible difference in geometrical acceptance between control samples (recorded 
with blocked beam) and signal samples (recorded with beam). The choice of control samples, 
selection and reconstruction procedures have been optimized to minimize these differences. 

Uncertainties on the muon halo background estimates are due to the limited size of the 
control samples, the uncertainty of their normalization caused by decays of beam kaons and 
pions upstream of the decay volume, and the time dependence of the geometrical acceptance. 

3.3.2 Kf^2 background in the Ke2 sample 

The K^2 background results mainly from muon misidentification {E/p > 0.95) due to 'catastrophic' 
bremsstrahlung in or immediately in front of the LKr calorimeter. It has been addressed by a 
dedicated measurement of the misidentification probability P^e based on samples of muons from 
Kfi2 decays traversing the Pb bar collected simultaneously with the K^(Ph) data samples. 
These samples are statistically independent from the signal i^^2 samples, as the Pb bar region 
is excluded from the geometrical acceptance, as discussed in Section [3.21 As noted in Section [21 
electrons passing through the Pb bar have a probability of only ~ 5 x 10"^ of being identified. 
Thus samples of misidentified muons traversing the Pb bar are electron-free for all practical 
purposes. However, the misidentification probability {Pj^g) for a muon traversing the Pb bar 
differs from the probability (P^e) for an unintercepted muon because of ionization energy loss 
(dominant at low momentum) and bremsstrahlung (dominant at high momentum). 
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Figure 3: (a) Misidentification probability for muons traversing the lead bar, Pj^^, for (£'/p)min = 
0.95 as a function of momentum: measurement (solid circles with error bars; the uncertainties 
are uncorrelated) and simulation (solid line), (b) Correction factors /pb = P^e/P^e evaluated 
by a simulation for the two specified values of {E/p)^ia. Dotted lines in both plots indicate 
the estimated systematic uncertainties arising from the simulation. The correlation of the latter 
uncertainties across the momentum values leads to a consistent shift of the MC band with respect 
to the measurements; the data/MC agreement validates the assigned systematic uncertainties. 



To evaluate the corresponding momentum-dependent correction factor /pb = Pfj,e/P^e 
measured probability Pj^^ for muons traversing the Pb bar, a dedicated MC simulation based on 
Geant4 (version 9.2) [13j has been developed to describe the propagation of muons downstream 
of the spectrometer involving all electromagnetic processes including muon bremsstrahlung ac- 
cording to Ref. [H]. The LKr calorimeter reconstruction has been optimized for showers ini- 
tiated by electrons and photons and starting near its front surface, whereas showers initiated 
by muon bremsstrahlung start throughout the detector volume. The corresponding systematic 
uncertainties on P^g and Pj^^ from the simulation due to energy calibration and cluster re- 
construction have been estimated to be 10% of their values. However, the uncertainty on the 
ratio /pb = Pi_ie/P^e is significantly smaller ((5/pb//pb = 2%) due to a cancellation of the main 
systematic effects. Moreover, /pb has low sensitivity to the model of muon bremsstrahlung 
cross-section used. The measurements of Pj^^ in momentum bins compared to the results of 
the MC simulation and the correction factors /pb obtained by the simulation, along with the 
estimated systematic uncertainties of the simulated values, are shown in Fig. [3l 

The K^2 background contamination due to muon misidentification has been computed using 
the geometrical acceptance evaluated by simulation, the measured P^^ and the correction for 
the effect of the Pb bar, /pb, evaluated from the simulation described above. The uncertainty 
on the background estimate comes from the limited size of the data sample used to measure Pj^^ 
and the uncertainty (5/pb of the MC correction /pb for ionization and bremsstrahlung in the Pb 
bar. Moreover, the positive correlation between the reconstructed M^^^^{e) and E/p, which are 
both computed using the reconstructed track momentum, leads to an apparent dependence of 
P^e on Af^jgg(e). The corresponding correction is based on the knowledge of the muon energy 
deposition spectrum in the vicinity of E/p = 1; its model dependence leads to an additional 
systematic uncertainty on the Kfj_2 background. The stability of the result with respect to 
variations in the lepton identification procedure was checked as discussed in Ref. [8]. 
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The K^2 decay also contributes to the background via /x^ — t- e^VeV^ decays in flight. This 
background has been evaluated with a MC simulation. Energetic forward secondary electrons 
compatible with K(,2 kinematics and topology are strongly suppressed by muon polarisation ef- 
fects and are further suppressed by radiative corrections to the muon decay [15] . The background 
from muon track association with accidental LKr clusters has been measured to be negligible. 

3.3.3 — )■ 6^1/7 background in the Kf,2 sample 

The definition of Rk includes IB but excludes structure-dependent (SD) radiation [1]. The 
SD'^ (positive photon helicity) component of the — ?■ 6^1/7 process peaks at high electron 
momentum in the rest frame (i?* ~ Mx/2) [12]; it is therefore kinematically similar to K(.2 
decay. It can contribute to the background if the photon escapes the acceptance of the LKr 
calorimeter. The background has been estimated by a MC simulation based on the measured 

—7- e^z^7 (SD"*") differential decay rate in the kinematic region E* > 200 MeV [7j. The main 
uncertainty on this estimate is due to the limited experimental precision on the decay rate. 

The — )• e^z^7 (SD~) decay with negative photon helicity peaking at E* « Mk/^ is 
kinematically incompatible with K(,2-, and the corresponding background is negligible. Similarly, 
the background from interference terms between the IB and SD processes is negligible. 

3.3.4 Other backgrounds in the Kf>2 sample 

The —7- TT^e^u and — t- tt^tt" decays produce a Ke2 signature in two cases: a) all vr*^ 
decay products are undetected and, for the latter decay mode, the tt^ is misidentified as e^; b) 
the only reconstructed particle is an electron (e^) from a Dalitz decay vr^, — )• e~^e~j. 

Due to the significant missing mass, these decays can only be kinematically compatible with 
Ke2 if the kaon is in the high-momentum tail of the beam distribution, or the detected kaon 
decay daughter particle suffers large multiple scattering. The systematic uncertainties on these 
minor backgrounds are due to the limited precision of the simulation of the non-Gaussian tails 
of multiple scattering; they have been estimated as 50% of the contributions themselves. 

The estimation of the — t- tt^tt^ background involves the pion misidentification proba- 
bility, which has been measured as a function of momentum from samples of — ?■ tt^tt^ and 

— )• TT^e^u decays (the latter collected during a special run). In particular, at high lepton 
momentum where — )• Tr^vr*^ contributes, the selection criterion 0.95 < E/p < 1.1 leads to a 
misidentification probability Pj^e = (0.41 it 0.02)%. 

Data samples collected with simultaneous and K~ beams, namely the K~{Ph) sample 
and a part of the K^(Ph) sample, are affected by backgrounds due to decays of opposite sign 
kaons, in a way similar to the muon halo control samples described in Section [3.3.11 Contribu- 
tions from the — )• n^i^u and — )• Tr^Tr*^ decays with subsequent vr^, — )• 7e"^e~ decays have 
been identified and subtracted using MC simulations. 

3.3.5 Summary of backgrounds 

The M^jgg(£) spectra of the selected K£2 candidates summed over all data samples are presented 
in Fig. SI The numbers of the selected K£2 candidates and backgrounds in data samples inte- 
grated over the lepton momentum are summarized in Table [H The total Ke2 sample consists 
of 145958 candidates with an estimated background of (10.95 it 0.27)%. The dependences of 
the backgrounds on the lepton momentum for the -R'^(noPb) and K^(Ph) data samples, which 
have the lowest and the highest background contaminations, respectively, are displayed in Fig. [5l 
The total Kf^2 sample collected with a pre-scaled trigger consists of 4.282 x 10^ candidates with 
a background due to the beam halo muons measured to be (0.50 ± 0.01)%. 
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Figure 4: (colour online). Distributions of reconstructed squared missing masses (a) Afmiss(^) 
and (b) M^-^^^[fi) of the candidates compared with the sums of normalized estimated sig- 
nal and background components. Beam halo contributions have been measured as discussed 
in Section 13.3. H the other contributions have been estimated with MC simulations involving 
the measured particle misidentification probabilities. The double peak structure of the -fC^2 
background in the Ke2 sample originates from the momentum dependence of the electron iden- 
tification condition. 

3.4 Other systematic effects 

A detailed discussion of the systematic effects not related to background subtraction is available 
in [8]. The main points are summarized below. 

3.4.1 Acceptance correction 

The typical geometric acceptance of the adopted K^2 selection is about 34% (42%) for the 
data samples collected with (without) the Pb bar. The acceptance correction A{Kfj_2)/A{Ke2) 
for each subsample, varying from 1.15 to 1.42 over the subsamples, has been evaluated by 
the MC simulation described in Section 13.11 Part of the correction, particularly at low lepton 
momentum, arises because of the wide radial distribution of electrons from Ke2 decays (not 
fully contained within the geometric acceptance of the DCH planes). Radiative effects result in 
a further loss of Ke2 acceptance by increasing the reconstructed M^-^^^{e) as quantified below. 
The radiative process — ?• e^vj from IB, which decreases the Ke2 acceptance by about 5% 
relative, is included in the simulation following \12\ I16j. The evaluation of the correction for the 
external bremsstrahlung, which accounts for about 6.5% relative Ke2 acceptance loss, requires a 
precise description of the material upstream of the spectrometer magnet jlOj : a Kevlar® window 
separating the vacuum decay volume from the spectrometer (0.3%Xo), helium in the 12.3 m long 
volume (0.4%Xo) and the two DCHs (0.8%Xo). 

The validity of the material description has been demonstrated by measuring the bremsstrah- 
lung rate for a sample of — )• 71*^6^ 1^(7) decays collected concurrently with the main data set. 
The fraction of events in this sample with a reconstructed bremsstrahlung photon of at least 
10 GeV energy, has been used to measure the material thickness. The quantity ranges 
from 0.3% to 2% for the momentum ranging from 20 GeV/c to 40 GeV/c, with contributions 
from external bremsstrahlung (~ 70%) and — )• ir^e^iyj (IB) decays (~ 30%). The relative 
uncertainty on the material thickness simulation has been conservatively estimated to be 1.1% by 
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Figure 5: The numbers of Ki2 candidates and main background contributions in lepton momen- 
tum bins in the i<r^(noPb) and i^^(Pb) data samples, which have the lowest and the highest 
background contaminations, respectively. The differences between the data samples in terms of 
the background composition are clearly visible. 
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Table 1: Background contaminations in the samples integrated over lepton momentum. The 
uncertainties on the background in the samples are negligible. 



Data sample 


i^+(noPb) 


A'+(Pb) 


ir-(noPb) 




Ke2 candidates 


59813 


63282 


10530 


12333 


Muon halo 


(1.11 ±0.09)% 


(1.51 ±0.10)% 


(4.61 ±0.18)% 


(7.86 ± 0.23)% 




(6.11 ±0.22)% 


(5.33 ± 0.19)% 


(5.76 ± 0.20)% 


(4.87 ±0.17)% 


K^2 (a* e decay) 


(0.26 ± 0.04)% 


(0.27 ± 0.04)% 


(0.31 ±0.09)% 


(0.19 ±0.07)% 


e^u-f (SD+) 


(1.07 ±0.05)% 


(4.01 ±0.18)% 


(1.25 ±0.06)% 


(3.95 ±0.17)% 




(0.05 ± 0.03)% 


(0.28 ±0.14)% 


(0.09 ± 0.05)% 


(0.37 ±0.17)% 




(0.05 ± 0.03)% 


(0.18 ±0.09)% 


(0.06 ± 0.03)% 


(0.18 ±0.09)% 


Opposite sign K 




(0.04 ±0.01)% 




(0.25 ± 0.03)% 


Total background 


(8.65 ±0.25)% 


(11.62 ±0.33)% 


(12.08 ±0.29)% 


(17.67 ± 0.39)% 


K^2 candidates /lO'' 


18.027 


18.433 


3.069 


3.288 


Muon halo 


0.39% 


0.44% 


0.77% 


1.22% 



comparing the values of measured from the data and from simulated — )• 7r^e^z^(7) samples. 
The quoted uncertainty includes contributions from the limited size of the control sample and 
the precision of the — )■ -K^e^v^ (IB) decay simulation according to [16]. It translates into 
a systematic uncertainty of 5Rk = 0.002 x 10~^ on the final result. No statistically significant 
variation of the material thickness over the data taking period (e.g. due to the helium purity 
variation) has been observed. 

An independent validation of the material simulation has been performed using special data 
samples taken on two earlier occasions with low-intensity, mono-energetic beams steered into 
the spectrometer. The rates of the radiated photons (produced only by external bremsstrahlung) 
agree between data and simulation within 1% precision. 

Other sources of systematic uncertainty are the limited knowledge of the beam spectrum, 
profile and divergence, and the simulation of the LKr response to soft radiative photons. A 
separate uncertainty has been assigned to account for the finite precision of the spectrometer 
alignment. 

3.4.2 Electron identification efficiency 

Samples of electrons and positrons selected kinematically from — t- i^^e^v and — t- -K^e^v 
decays have been used to calibrate the energy response of each LKr cell for each of the nine data 
taking periods defined for this purpose, as well as to measure the electron identification efficiency 
/e as a function of LKr cell, track momentum and time. The overall inefficiency averaged over 
the Ke2 sample is 1 — /e = (0.72 ± 0.05)%, where the quoted uncertainty is mainly systematic. 
The corresponding correction to the final result is ARk = (±0.018 ± 0.001) x 10~^. 

The inefficiency has been found to be highly stable in time: averaged over the momentum 
ranges p < 25 GeV/c {p > 25 GeV/c), it varies during the data taking from 0.43% to 0.44% 
(from 0.79% to 0.81%). Note that the electron identification criteria are different in the two 
momentum ranges, as explained in the last paragraph of Section [3.2i 

3.4.3 Trigger and readout efficiencies 

The (1.4 ±0.1)% inefficiency of the Qi trigger condition, as well as the much smaller inefficiency 
of the 1-track trigger condition, nearly cancel between the Ke2 and K^2 samples due to their 
geometrical uniformity and long-term stability. The residual systematic bias is negligible. The 
inefficiency of the Slki condition, which enters only into the Ke2 trigger chain, has been mea- 
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Figure 6: Stability of the Rk measurement versus lepton momentum and for independent data 
samples. Statistical and total errors of individual points are indicated with horizontal dashes. 
The result of the fit to the 40 subsamples, with the combined statistical and systematic error, 
is shown by the horizontal lines and bands. 



sured with a sample of — t- events. It is found to be significant only in the lowest bin 

of lepton momentum (close to the 10 GeV trigger energy threshold). It varies in that bin over 
time from 0.15% to 0.5%, and is localized geometrically. The correction to the final result for 
the -BLKr trigger inefficiency amounts to ARk = +0.001 x 10~^, with a negligible uncertainty. 

Energetic photons not reconstructed in the LKr may initiate showers by interacting with 
the DCHs or the beam pipe material, causing the DCH hit multiplicities to exceed the limits 
allowed by the 1-track trigger condition. The dominant background with a lost photon is from 
the — )• e^i^7 (SD"*") decay. The corresponding 1-track trigger inefficiency for the background 
— )• e^uj (SD"*") events in the Ke2 sample has been evaluated by a MC simulation. The 
simulation has been validated by comparison with a data sample of — ?• -K^e^v events with 
two lost photons. The inefficiency has been found to vary from 0.06 to 0.21 depending on the 
electron momentum, and its relative systematic uncertainty has been estimated to be 30%. The 
correction to Rk for this effect is ARk = (+0.003 ± 0.001) x 10"^ 

The misalignment of the LKr signal timing leads to a small global LKr readout inefficiency 
affecting the Ke2 reconstruction only. This inefficiency has been measured to be 1 — /LKr = 
(0.20 ± 0.03)% and found to be stable in time, using an independent readout system with 
larger granularity developed to monitor the trigger chain of the NA48 experiment |10] . The 
corresponding correction to the result is ARk = (+0.005 ± 0.001) x 10~^. 



4 Averaging over the data samples 

A fit to the 40 measurements of Rk in individual subsamples (consisting of the 4 data 
samples, each subdivided into 10 lepton momentum bins) has been performed. The following 
correlations of the systematic uncertainties over subsamples have been taken into account. 

• All systematic uncertainties related to the acceptance correction and electron identification 
efficiency are conservatively considered to be fully correlated among data samples and 
lepton momentum bins. 
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Table 2: Summary of the uncertainties on Rk- 



Source 


5Rk X 10^ 


Statistical 


0.007 


K^2 background 


0.004 


— )• e^uj (SD"*") background 


0.002 


—7- TT^e^u, — >• TT^vr^ backgrounds 


0.003 


Muon halo background 


0.002 


Spectrometer material composition 


0.002 


Acceptance correction 


0.002 


Spectrometer alignment 


0.001 


Electron identification inefficiency 


0.001 


1-track trigger inefficiency 


0.001 


LKr readout inefficiency 


0.001 


Total systematic 


0.007 


Total 


0.010 



• The systematic uncertainties on the decay rates of all background kaon decay modes, com- 
ing from other measurements [iTj , are fully correlated among data samples and lepton mo- 
mentum bins. The most significant of these uncertainties is that of the — )• e^v^ (SD"*") 
decay rate. 

• The systematic uncertainties on the muon and pion misidentification probabilities have 
components which are fully correlated among the data samples (the statistical errors of 
independent P^^ and P-,^e measurements in a given lepton momentum bin) as well as 
components which are fully correlated among both data samples and lepton momentum 
bins (due to the model-dependent correction described in Section I3.3.2p . 

• The systematic uncertainties on the muon halo background normalization, which is relevant 
only for high lepton momentum bins, are fully correlated among the data taking periods. 
Uncertainties due to the residual differences of geometrical acceptances among data and 
muon halo control samples are correlated among lepton momentum bins. 

• The systematic uncertainties due to the 1-track trigger inefficiency and LKr readout inef- 
ficiency are correlated among the data samples. 

• The systematic uncertainties on the — )• n^e^u and — )• tt^tt^ backgrounds due to 
imperfect simulation of multiple scattering at large angles and of the kaon beam spectrum, 
which affect the highest lepton momentum bin only, are correlated among the data samples. 

The resulting correlation coefficients among the total uncertainties for the 40 subsamples range 
from 0.002 to 0.40. The fit result is 

Rk = (2.488 ± 0.007stat. ± 0.007syst.) x 10"^ = (2.488 ± 0.010) x 10'^, 

with x^/ndf = 47/39 (probability 18%). The stability of the measurements for bins of different 
lepton momentum (averaged over the four data samples) and for the different data samples 
(averaged over the ten momentum bins) is shown in Fig. [H The contributions to the uncertainty 
of the result are summarised in Table [2j 
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Summary 

The most precise measurement of Rk = T(Ke2)/^{K^2) to date has been performed from a 
sample of 145958 Ke2 candidates with an estimated background of (10.95 ± 0.27)% cohected by 
the NA62 experiment in 2007-2008. The result Rk = (2.488 ± 0.010) x 10"^ is consistent with 
the earlier measurements and with the SM expectation. The experimental uncertainty on Rk is 
still an order of magnitude larger than the uncertainty on the SM prediction, which motivates 
further measurements at improved precision. 
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